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ABSTRACT 

We propose that bound, young massive stellar clusters form from dense clouds that have escape speeds 
greater than the sound speed in photo-ionized gas. In these clumps, radiative feedback in the form of gas 
ionization is bottled up, enabling star formation to proceed to sufficiently high efficiency so that the result- 
ing star cluster remains bound even after gas removal. We estimate the observable properties of the massive 
proto-clusters (MPCs) for existing Galactic plane surveys and suggest how they may be sought in recent and 
upcoming extragalactic observations. These surveys will potentially provide a significant sample of MPC can- 
didates that will allow us to better understand extreme star-formation and massive cluster formation in the Local 
Universe. 

Subject headings: galaxies: stars: massive, - star formation: cluster formation, - HII regions; - ISM: bubbles 



1. INTRODUCTION 

The formation of bound star clusters has become a topic of 
renewed interest. The Milky Way contains about 150 globular 
clusters (GCs) with masses from 10 4 M Q to over 10 6 Mq and 
tens of tho usands of open clusters conta ining from 100 to over 
10 4 stars (Porte gies Zwart et al.||2010| . While no GCs have 
formed in the Milky Way within the last 5 Gyr, bound clusters 
that survive for more than hundreds of crossing times continue 
to form. 

Infrared observations over the last two decades have shown 
that molecular clouds tend to produce stars in highe r surface 
densities (> 3 stars pc~ 2 ) than the field population ( Lada & 



|L,ada||2003) . |Bressert et~aT] ( |2010) > showed that stars within 
500 pc of the Sun form in a smooth continuous distribution 
and only a minority will dynamically evolve to form bound 
low-mass stellar clusters (10 2 to 10 3 M Q ). The vast majority 
of these young clusters are transient groups that are bound pri- 
marily by the gas in their environment. Thus, while most stars 
may form in groups, gravitationally bound clusters which re- 
main bound for many crossing-times following dispersal of 
their natal clump are rare and contain less than 10% of the 
Galactic stellar population. Despite the small number of stars 
that fo rm in the bound young m assive clusters (YMCs; > 10 4 
Mq ; |Portegies Zwart et al.|2010 1, they are important as they 
shed light on extreme star-formation in the Local Universe 
and provide insight on how GCs may have formed in the high- 
redshift universe and in the distant past of the Milky Way 
(Elmegreen & Efremov 1997). 

With new Galactic pla ne surveys, e.g., the Herschel Hi- 
GAL survey ( |Molinari et al.||2010) , the AP EX Telescope 
Large Area Survey of the Galaxy (ATLASGAL;|Schuller et al.| 



2009), the Bolocam Galactic Plane Survey (BGPS; Aguirre 
et al.|2011), the H 2 Q Southern Galactic Plane Survey (HOPS; 
Walsh et al.||2011|), and the Millimeter Astronomy Legacy 



Team 90 GHz Survey (MALT90; |Foster et al.|20~TT) , we are on 
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the cusp of better understanding how massive clusters form. 
The question is, how do we find the massive proto-clusters 
(MPCs) that will form these YMCs? We investigate how the 
YMCs may form and provide a simple model with observa- 
tional properties that can be used to identify MPC candidates. 
To identify more extreme MPCs in nearby galaxies (e.g., the 
Antennae Galaxies), we need capable telescopes like the Ata- 
cama Large Millimeter/sub-millimeter Array (ALMA). 

YMCs may predominantly form from MPCs having gravi- 
tational escape speeds greater than the sound speed in photo- 
ionized gas. When this condition is met, ionization cannot 
disrupt the entire MPC. Stars can continue to form from the 
remaining neutral gas and star formation efficiency (SFE) in- 
creases to 30% and higher. The remaining mass in the stellar 
population nullifies the effects of gas expulsion and the cluster 
will remain bound. If the absolute value of the gravitational 
potential energy is greater than the expected thermal energy of 
the plasma in a massive gas clump, and supemovae have not 
yet occurred, then we consider the object to be an MPC can- 
didate. We compare the effective photo-ionized sound speed 
of the plasma (c n ) to the escape velocity (v esc ) of the clump, 
which we denote as SI = v esc /c n . A gas clump that has S7 > 1, 
implying v esc > c u , is an MPC candidate while SI < 1 are not, 
since gas can be dispersed by the appearance of the first OB 
stars. 

We describe the simple model of massive stellar cluster for- 
mation using the ft parameter in § [2] We make predictions 
on the MPCs observational properties for Galactic plane sur- 
veys and ALMA (extragalactic) in §[3] §|4] discusses the im- 
plications of the model and predictions with a summary of the 
results. 

2. YOUNG MASSIVE CLUSTER FORMATION 

2.1. Initial Conditions 

What are the initial conditions that lead to bound stellar 
cluster formation, e.g., YMCs like NGC 3603 and R136? We 
must consider the differences between low-mass and high- 
mass star forming clouds. In the solar neighborhood the local 
SFE, defined as the final mass in stars formed in a cloud di- 
vided by the initial mass of gas, of th e low-mass star fo rming 
regions is reported to be 5 % or less (|Evans et al.|200"9]). Nu- 
merical experiments (e.g., |Lada et al.|1984||Geyer & Burkert 
|2001 1 [Goodwin & Bastian 2006) demonstrate that the forma- 
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tion of gravitationally bound clusters requires that the local 
SFE ^30%. Observations show similar results of high er SFE 
for massive star forming regions(e.g., |Adams||2000l |Nurn 



|berger et al.|2002| l and hint that some of these systems could 
evolve into open clusters. The initial conditions that enable 
the formation of YMCs may be responsible for the high SFE 
observed in these bound stellar clusters. From the gas clump 
phase to the stellar clusters the gas is somehow efficiently con- 
verted to stars. Assuming that the SFE of 30% is critical for 
forming a YMC, then the MPC mass should be greater than 
3xlO 4 M . 

Star formation in a clump ends when its gas is dispersed by 
the feedback energy injected from its newborn stars. Gas can 
be dispersed when the outward pressure generated by feed- 
back exceeds the inward pressure of the overburden of gas in 
the cluster gravitational potential well (Bally 2011). The ra- 
tio of the local escape velocity divided by the sound speed for 
photo-ionized gas, f2 = v esc /c n , plays a crucial role in deter- 
mining if a clump forms a bound cluster or not. Feedback 
from accre ting low-m ass protostars is dominated by bipolar 
outflows (Ballyi [20TT) . As stars reach several solar masses, 
their non-ionizing near-UV radiation photo-heats surrounding 
cloud surfaces, raising the sound speed to ^ 1 to 5 km s~ . 
Stars greater than 10 M Q (early-B and O stars) ionize their 
surroundings, raising the sound speed up to ~ 10 km s -1 . Al- 
though protostellar outflows, stellar winds, and radiation pres- 
sure of ever-increasing strength can also raise the effective 
sound speed by generating internal motions, only the effects 
of ionizing radiation will be considered here. See § 2.2 for 
details. 

If f2 < 1, the gas can be dispersed in a few crossing times 
from the star-forming clump, bringing star formation to a halt 
with a low stellar density. On the other hand, when il > 1 
the gas will remain bound and can continue to form new stars 
or accrete onto existing ones until further increases in stellar 
luminosity or mechanical energy injection raise the effective 
sound speed to a value greater than the escape speed. As the 
stellar mass increases , energy released by t he forming em- 
bedded cluster grows ( Miesch & Bally|1994[ l. In such a clus- 
ter, photo-ionized plasma will remain gravitationally bound 
by the cluster potential and recombinations in the ionized 
medium will tend to shield denser neutral clumps, allowing 
star formation to proceed to high efficiency. This model works 
under the assumption that the entire gas clump is instantly 
ionized, which is a worst case scenario for a clump to remain 
bound. This means that even if there is a large number of 
OB stars present in a clump, its gas will unlikely be fully ion- 
ized and our model would still hold. We discuss the details of 
ionization and its effect on pressure balance and ongoing star 
formation below. 

2.2. Ionizing Feedback, Pressure, and Star Formation 

Feedback, which plays a major role in the self-regulation 
of star formation can be subdivided into two forms: mechan- 
ical feedback consisting of protostellar outflow, stellar winds, 
and supernova explosions and radiative feedback which can 
be subdivided into radiation pressure, non-ionizing FUV heat- 
ing, and ionizing EUV heating. In massive cluster forming 
environments, prior to the explosion of the most massive star, 
stellar winds will dominate mechanical feedback and all three 
radiative mechanism can be important (Bally 201 1 ). 

Consider a fiducial reference clump of M gas = 2 x 10 4 M Q 
and Madia, =~ 1 x 1O 4 M (SFE = 30%), a total luminos- 



ity Lf, Lyman continuum photon luminosity Qf, and a total 
stellar wind mass loss rate Mf at Vf (wind terminal veloc- 



ity) which corresponds to 50 07 stars (30 M ) located at the 
center of the cluster. The assumed quantities for the fiducial 
cluster ar e shown in the denominators of the eq uations be- 
low (see |Martins et al.|[2003) |Donati et al.|[2 002, and refer- 
ences therein). Through 10 4 Monte Carlo iterations the ratios 
of L, Q, and M between the stellar population of the fidu- 
cial cluster and Salpter-like cluster (8 M Q < M slar < 1 10 M Q ) 
are shown to be L f /L Salpeter = 14.0, Qf/Q Sa i peter = 2. 2, and 
M f /M Salpe 



■ 8.9. These calculations were done using 



Crowther 



Mur- 



(2000) for M 

2 
II 



The 
2.13 pc 



ray & Rahman| ( |2010| ) for Q and|i_ 
gravitational radius for this cluster tq = 2GM/cj 
and as ~ 2.6 x 10~ 13 cm 3 s _1 is the case-B recombination rate 
for H at a temperature of 10 4 K. 

Below, we consider the radiation pressure, stellar wind ram 
pressure, and the internal pressure of a uniform density H II 
region. These pressures are evaluated at D = 2 pc from the 
center of the cluster, a distance close to the gravitational ra- 
dius for the reference cluster. 



1 . Stellar wind: 
P w = p (r)V 2 = MV /AttD 1 
= 1.32x 10" 10 dynes cm -2 



Mf 

5xlO- 6 M yH 



2x 10 3 km s" 



1 pc 



(1) 



2. Radiation pressure on an optically thick surface: 
P md = L/4ircD 2 

= 2.68 x 10~ 9 dynes cm" 



-2 






\Df] 




_ W L Q _ 




2 pc 
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3. Thermal pressure of a uniform density H II region with 
Stromgren radius equal to D = 2 pc: 



Pn u = pLm H cl(3Q/4ira B y/ 2 D- 3 / 2 



■■ 5.40 x 10 dynes cm 



-2 


Qf 


1/2 


El 




10 51 photons s 1 




2 pc 



-3/2 

(3) 



These pressures are listed in order of increasing significance 
for the fiducial reference cluster. It is important to note that 
the equations express the feedback pressures as power laws of 
distance from the center of the fiducial cluster. The pressure 
from Eq. 3 tapers off the slowest among the three pressures, 
which bolsters the dominance of photo-ioniza tion at large dis- 
tances . A simil ar result was fo und in lKrumholz & Matznerl 
(2009), but Murray et al. (20101) conclude that radiation pres- 



sure is the dominant mechanism in dispersing gas under dif- 
ferent scaling assumptions. If radiation does play a role in dis- 
persing the gas it could mass-overload the shell around the H 
II region and force some gas to escape, which may imply that 
a YMCs SFE is unlikely to ever be 100%. However, the ion- 
ized gas pressure will not unbind the system and could create 
large optically thin bubbles, allowing much of the radiation 
pressure to escape without interacting with any gas particles. 
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3. PREDICTING OBSERVED PROTO-CLUSTER PROPERTIES 

3.1. Proto-cluster Geometries 

What are the expected physical geometries of MPCs? Us- 
ing both the YMC progenitor model and observations of 
YMCs, we derive the maximum sizes and mass of the MPCs. 
In the bound H II gas clump model, the key requirement is 
v esc > c n to form a YMC. By setting the potential and kinetic 
variables to —GMmr~ x and \mc\, we can solve for r. We de- 
note this radius as tq, which is a function of the clump mass, 
r n(M c i wnp ) = 2GM c i ump c~^ . For MPC candidates with masses 
of 3 x 10 4 - 3 x 10 6 Mq, rn ranges between 2.6 and 258.1 pc 
(see Table 1 for details). The rn values for > 10 6 Mq clumps 
are large and are very unlikely to form a > 10 6 Mq gravita- 
tionally bound star cluster as the free-fall time of the system 
could exceed when supernovae could begin and disrupt the 
system. Hence, an upper size limit for these objects is needed. 

There is a well-measured constraint on YMCs that we can 
apply to predict the upper limit radii for the MPCs. Re- 
cent high resolution imagining and spectral studies of YMCs 
like NGC 3603, Arches, Westerlund 1, and R136 have shown 
these systems t o be in or close to virial equ ilibrium at ages 
of ~ 1 Myr dRochau et al.|[20T0l |Clarkson et al. | [20TT] |Cot"| 
|taar et al |[2012t |Hen ault-Bru netet al.||2012| >. This impE 



that the YMCs must have gone through at least one full 
crossing time before their pre sently observed age. U s ing th e 
crossing time equation from Portegies Zwart et al. ( 20101 , 
tcmss = (GMr~?) -1 ' 2 , and fixing the crossing time to 1 Myr we 
can solve for r n > = (GMf^^) 1 / 3 , which is only dependent on 
the mass. For the same clump mass range mentioned for 
above, r v( > spans from 5.1 to 23.8 pc (see Table 1 for details). 
Note that the rn and r v( > values are the upper limit radii for the 
MPCs. This is an important aspect to keep in mind as there 
is no evidence for YMC s to have a pro portionality between 
mass and radius (Larsen 2004; Ba stian et al.|2012| ). rn and 
r VJ > have similar radii between 10 4 M e and 10 Mq intersect- 
ing at 8.4 x 1O 4 M , but at > 10 6 Mq r n > r vir . Extragalactic 
predictions regarding more massive MPCs (3 x 10 7 - 3 x 10 9 
Mq ) is provided in Table 1 . 

Figure 1 shows the upper limit radii relative to clump 
masses for rn and r,,,>. We include a shaded region mark- 
ing where MPC candidat es reside. Infrared dark clou ds and 
clump s reported in Rath borne et al.| ( |2006[ l and [Walsh et al.| 
(201 1 ) are included to show where they lay on the plot relative 



to rn and r V!r . Five g as clumps fit within the MPC criterion: 
GO.253+0.016 ( [Longmore et al.|2012b, an MPC candidate in 
the Antennae galaxy (Herrera et al. 2012), and three Galactic 
clumps reported in Ginsburg et al. (2012) which are discussed 
in further detail in their paper. 

We took a collection of G alactic YMCs summarized in 
Portegies Zwart et al. ( |2010| l and estimated their progenitor 



masses by multiplying their current stellar mass by a factor 
of 3 to account for SFE^ 30%. Eight of the of 12 estimated 
YMC progenitors are within the virial and f2 radii limits. The 
other four are found not only close to the upper limit of the rn 
line, but appear to follow the line. If we assume that the SFE 
> 30% then the estimated YMC progenitors will move above 
rn line. 

3.2. Prediction and Observations 

From the derived physical properties of the MPCs dis- 
cussed, we predict the MPCs integrated fluxes in the HiGAL, 
ATLASGAL, and BGPS surveys in Table 1 using mass to flux 



conversions discussed in Kauffmann et al. (2008 ) assuming T 
= 20 K. If T = 40 K, then the flux will increase by a fac- 
tor of ~ 2. Between the wavelengths of 500 pm and 1 100 
the emission from the bright MPCs and the interstellar 
medium are optically thin throughout the Galactic plane. We 
predict that the surveys should be sensitive to all of the MPCs 
in the Milky Way. Tackenberg et al. ( 2012[ ) independently 
came to similar conclusions for ATLASGAL. Moreover, the 
sources for the mass clumps > 3 x 10 4 Mq are resolvable up 
to 20 kpc. We calcu lated the clump central pressures from 
Johnst one et al.J ( [2000| . These values are provided in Table 
1. Note that HiGAL has advantages over ATLASGAL and 
BGPS, since it has no spatial filtering and has the highest sen- 
sitivity amongst these surveys. 



Ginsburg et al. (2012 1 used the predicted properties to de- 
velop a selection criteria from the BGPS catalog and discov- 
ered three MPC candidates. The masses of the objects can be 
determined from the dust emission (Mj„ s ,), but comparing the 
dust masses to virial mass, M VJ >, is important in determining 
whether a MPC candidate is gravitationally bound or not. Fol- 
lowing similar treatment from Longmore et al. ( 2012| l one can 
use AV from molecular line tracers (e.g., HCO+ and N2H+) 
to approximate a clump's virial mass, M,,,> oc r c i ump AV 2 . If 
Mdusi ~ M v i r , then we can assume that the clump is consistent 
with being gravitationally bound. Estimating the mass of a 
clump from continuum emission requires a known distance to 
the objects using Vlsr, which line surveys can provide (e.g., 
Walsh eTal.|201 1| Schlingman et al.|201 l| l. The combination 



of these line and continuum Galactic plane surveys will allow 
us to obtain a near complete census of MPCs in the Milky 
Way. 

In Table 1, we show our predictions for the most massive 
MPCs (> 10 7 Mq) in the Local Universe regarding ALMAs 
observing capabilities. The given fluxes are calculated for the 
MPCs at a distance of 40 Mpc and 6 ~ 0".5. These objects, 
if near the upper limit of r V! >, could be resolvable at these dis- 
tances. 

4. DISCUSSION AND SUMMARY 

We have discussed the possible conditions necessary for 
YMC formation and how to identify their progenitors, MPCs, 
in their primordial state regarding their masses, radii, and flux 
brightness. The key to identifying whether a massive gas 
clump can form a YMC is the balance between the gravi- 
tational potential of the gas clump and the gas kinematics. 
We characterize this balance as the ratio between v esc at a 
given radius and the sound speed of the photo-ionized gas, 
Cu ~ 10km s~ . If fl = v esc /c u > 1 (equivalent to, v esc > c n ) 
for a clump of gas (> 3 x 1O 4 M ) then the system will op- 
timally convert the clump to stellar mass and likely form a 
YMC. We classify such clumps as MPC candidates. If fi < 1 
then the system does not have a deep enough potential well to 
keep the photo-ionized gas bound, which will lead to rapid gas 
dispersal and low star formation efficiency. The end product 
will be a low mass cluster or group that will feed its stars to 
the field star population over a short time scale. It may be pos- 
sible that some of these MPC candidates will not form YMCs 
due to the cruel cradle effect, wher e the forming cluster is 
disrupted by nearby massive GMCs ( jKruijssen et al.|201 l| l. 

With the HiGAL, ATLASGAL, and BGPS Galactic plane 
surveys we should be sensitive to the MPCs throughout the 
Milky Way. The combination of these surveys will provide us 
a near complete sample of the MPCs in the Milky Way and 
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FIG. 1. — The mass-radius parameter-space for clumps partitioned by radii for rn (solid blue) and r w > (solid black). MPC candidates are defined with the 
following properties (green shaded region): a minimum mass of 3 X 10 4 Mq, r > rn for M c \ umB < 8,4 X 10 4 Mq, and r > r H > forMcj„„,„ > 8,4 X 10 4 M p. Clump 
masses and sizes are plotted on top from three different data catalogs: IRDCs (Rathborne et al. 2006 1, HOPS clumps I Walsh et al. 201 1 1, and YMCs (Portegies 
| Zwart et al.|2010 ). The YMCs are converted to their possible clump progenitors by assuming that Shb is ~30%, which boosts the mass of the systems by a tactor 
ot 10/3. The scaled YMC progenitors happen to lie near the critical rn line with out any tweaking of parameters. Two published sources that have radii less than 
both their respective and r,,;, are GO. 253+0. 016 (L12; Longmore et al. 2012) and an extragalactic massive proto-cluster candidate reported in Herrera et al. 
(H12; 2012). The MPC candidates reported in Ginsburg et al. (Gt>; 2012) are shown as squares. 



TABLE 1 

Predicted Proto-cluster Properties 







Galactic 


d < 20 kpc 










Mass 


logPcenA 




' vir 


&20 kpc 


-F500 


^850 
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Mq 


K7m 2 


pc 


pc 




Jy 


Jy 


Jy 


3 x 10 4 


13.15 


2.6 


5.1 


53.62 


30.78 


4.17 


1.81 


3 x 10 5 


13.81 


25.8 


11.0 


226.89 


17.19 


2.33 


1.01 


3 x 10 6 


14.47 


258.1 


23.8 


490.91 


36.72 


4.98 


2.16 






Extragalactic 


d < 40 Mpc 


8 ALMA = 0"-5 








Mass 


logPcenA 




r v ir 


#40 Mpc 


^450 


^850 


-F1200 


M Q 


K7m 2 


pc 


pc 




mJy 


mJy 


mJy 


3 X 10 7 


15.13 




51.2 


0.53 


77.7 


7.8 


2.3 


3 x 10 s 


15.80 




110.3 


1.14 


167.5 


16.7 


5.1 


3x 10 9 


16.47 




237.7 


2.45 


360.6 


36.6 


11.0 



Top half: The predicted physical properties of the YMC progenitors, which are observable throughout the Milky Way (<=20 Kpc) with the Hi-GAL, ATLASGAL and BOLOCAM 
Galactic plane surveys. Bottom half: the predicted physical properties of the YMC progenitors, which are observable up to 40 Mpc with ALMA assuming an angular resolution of 0.5". 
The values for are not shown for the extragalactic sources since r Vi > <C in all cases and hence such objects are immediately considered as PMC candidates. All fluxes are derived 
using r^i for the < 8.4 X 10 4 Mq and r v; > for > 8.4 X 10 4 Mq . Note that the fluxes for the extragalactic clumps are calculated using the rest wavelengths. 
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pave a path for future high-resolution studies. Furthermore, 
with ALMA's full potential we should be sensitive to MPCs 
over > 10 7 Mq in extragalactic systems within 40 Mpc of 
the Milky Way. This would help us better connect the Galac- 
tic and extragalactic MPCs to better understand the formation 
and evolution of YMCs, some of which could produce long 



lived "young globular clusters" (see Portegies Zwart et al. 
2010 and references therein). It is important to note that r vlr is 



the upper limit to these very massive proto-clusters and such 
objects will likely not be close to such scales. Observing these 
extragalactic MPC candidates in the nearby galaxies will help 
constrain the upper limit radii. 

This work was supported by NSF grant AST0407356. The 
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